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Abstract. We discuss the polarization signature of primordial gravitational waves 
imprinted in cosmic microwave background (CMB) anisotropics. The high-energy 
physics motivated by superstring theory or M-theory generically yields parity violating 
terms, which may produce a circularly polarized gravitational wave background (GWB) 
during inflation. In contrast to the standard prediction of inflation with un-polarized 
GWB, circularly polarized GWB generates non- vanishing TB and EB-mode power 
spectra of CMB anisotropics. We evaluate the TB and EB-mode power spectra taking 
into account the secondary anisotropics induced by the reionization and investigate 
the dependence of cosmological parameters. We then discuss current constraints on 
the circularly polarized GWB from large angular scales (£ < 16) of the three year 
WMAP data. Prospects for future CMB experiments are also investigated based on a 
Monte Carlo analysis of parameter estimation, showing that the circular polarization 
degree, e, which is the asymmetry of the tensor power spectra between right- and 
left-handed modes normalized by the total amplitude, can be measured down to 
|e| >0.35(r/0.05)-°' 6 . 
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1. Introduction 

The gravitational wave background (GWB) originating from inflation provides a direct 
probe of inflation against which we can test inflationary models. In particular, the energy 
scale of inflation can be estimated from the amplitude of the GWB. Combining with 
observations of scalar-type fluctuations, the detection of the GWB directly constrains 
the inflaton potential [JJ. Currently, there is no rigorous constraint on the amplitude 
of the GWB characterized by the tensor-to-scalar ratio, r = Aj^(A;o)/A^(A;ojjl In the 
standard scenario of slow-roll inflation, the GWB is expected to have a nearly scale- 
invariant spectrum, suggesting that the GWB would be detectable in a wide range 
of wavelengths or frequencies. For a large-scale experiment, polarization anisotropies 
of the cosmic microwave background (CMB) would be a powerful tool to search for 
primordial tensor fluctuations. Indeed, post-Planck experiments such as SPIDER^ and 
CMBpol (or Inflation Probe) in the Beyond Einstein program of NASA0 is dedicated 
to measuring the B-mode polarization anisotropies originating from the inflationary 
GWB, with expected precision level r ~ 10~ 3 [21 [3]. On the other hand, a direct 
measurement of the stochastic GWB might be possible for a small-scale experiment, 
especially using space-based laser interferometers [U El El [7] . Proposed missions such 
as the Big-Bang Observer (BBO) [8] and the deci-hertz gravitational-wave observatory 
(DECIGO) [HI [10] indeed aim at detecting the primordial gravitational waves at the 
frequencies / ~ 0.1 — 1Hz. Notice that the observational frequencies (or wavelengths) of 
the space interferometers are greatly different from those of the CMB experiments by 16 
orders of magnitude. This implies that the combination of both experiments provides a 
stringent constraint on the dynamics of inflation. 

Meanwhile, motivated by high-energy physics, there are numerous discussions on 
the corrections to the prediction of standard inflationary models. For example, some 
inflationary models contain parity-violating interaction terms, as generic predictions of 
superstring theory/ M-theory [TTl [T2] . Among these, the Chern-Simons term, which is 
a higher-order curvature term coupled to the scalar field, appears through the Green- 
Schwarz mechanism and the cosmological implications of the Chern-Simons terms has 
been extensively discussed [131 HH EE ESI El EES EES ED] . Such parity- violating terms 
directly affect the tensor-type perturbations [2T1 122] and the polarization modes of the 
resulting tensor fluctuations becomes asymmetric, leading to a circularly polarized GWB 
[TBI dU [15l ED]- In this respect, the detection of a circularly polarized GWB would be a 
direct signature of the cosmological parity-violation and it might also imply that there 
is a fundamental theory of particle physics beyond the standard model Qfl. Since the 
sensitivity of the forthcoming CMB experiments will be improved significantly, it is 

| In this paper, we adopt the conventional value fco = 0.002Mpc -1 as the pivot scale. 
§ http: //www. astro . caltech. edu/ ^lgg/spider_f ront . htm| 
|| http : //universe . nasa . gov/program/ inflation . html 

% There exists another source to generate a circular polarized GWB in the early universe, i.e., primordial 
helical turbulence produced during a first-order phase transition |231 [24] . However, the wavelength of 
the produced GWB is much smaller than the scale of CMB observations. 
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timely to explore the possibilities of measuring the signature of new interaction terms 
through the CMB anisotropies. 

In this paper, we discuss in some detail, the observational possibilities of using CMB 
anisotropy measurements to probe the additional signature imprinted in the primordial 
gravitational waves. According to p3], a circularly polarized GWB produces a non- 
trivial correlation between the temperature and the polarization anisotropies. As a 
result, the cross power spectra between temperature and B-mode polarization becomes 
non- vanishing. They calculated the TB-mode spectrum in an idealistic situation with a 
large tensor-to-scalar ratio, neglecting the secondary anisotropies. In the present paper, 
extending their analysis, we quantitatively evaluate the TB mode spectrum taking into 
account the effect of secondary anisotropies. Also, we calculate another non-vanishing 
spectrum, the EB-mode spectrum. Based on the three year data of Wilkinson Microwave 
Anisotropy Probe (WMAP) [251 EE], we discuss the current constraint on the degree 
of polarization of the primordial GWB. Further, we address future prospects for the 
PLANCK satellite or cosmic-variance limited experiments and estimate the extent to 
which the degree of polarization is constrained from future observations. 

This paper is organized as follows. In £j2j we briefly describe a mechanism to 
generate circular polarization of the GWB through the gravitational Chern-Simons 
term. In §|3l the influence of the polarized GWB on the CMB anisotropies is discussed. 
Based on this, the CMB power spectra originating from circularly polarized GWB are 
calculated in §H The current constraints and future prospects are discussed in §|5j 
Finally, £0is devoted to discussion and conclusions. 

2. Polarized gravitational waves from gravitational Chern-Simons term 

In this section, we briefly review a mechanism to generate a circularly polarized GWB by 
the gravitational Chern-Simons (gCS) term [13]. In superstring theory/M-theory, there 
exist scalar fields coupled with anti-symmetric tensor F A F = e 01 ^ F a ^F^ and/or R A 
r = e al3l& R a p y R-yS^u, where F^ v is the field strength of the electromagnetic field, R a/ 3 7 s 
is the Riemann tensor, and e a/3lS is a totally antisymmetric Levi-Civita tensor density 
[TT| [12]. These two terms are referred to as the electromagnetic and the gravitational 
Chern-Simons term, respectively. The presence of such parity-violating terms plays 
an important role for several cosmological issues such as structure formation involving 
axions [H] and leptogenesis or baryogenesis in the early universe [161 EL71 CHI HH] • m a 
homogeneous and isotropic background spacetime, the electromagnetic Chern-Simons 
term affects neither the evolution of the background spacetime nor the evolution of 
fluctuations as long as only linear perturbation is considered [27]. Therefore, we consider 
the gCS term only: 



where m v \ is the Planck mass, and the function f(<j)) represents a generic coupling 
to the scalar field 0. In some cases, the scalar field is identified with the inflaton 




(1) 
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field during inflation. As long as the inflaton field is homogeneous and constant in 
time, equation ([1]) is just a surface term, and it does not contribute at all to classical 
gravitational dynamics. Thus, after the end of inflation, we expect that the classical 
evolution without gCS term is recovered and no anomalous parity violation appears. 
Moreover, the gCS term also does not affect the evolution of the background and scalar 
perturbations in the linear regime [211 122] • Thus, if we ignore the vector perturbation, 
the influence of the gCS term only appears in the evolution of tensor perturbations. 

Let us linearize the Einstein-Hilbert action in the presence of the gCS term. 
Assuming a flat Friedmann-Robertson- Walker cosmology, the corresponding metric 
neglecting the scalar perturbation takes the following form: 

ds 2 = a 2 (r])[-dr] 2 + ((% + h ij )dx i dx j ] , (2) 

with hij being a transverse and traceless tensor, i.e., d^hij = h l i = 0. Expanding the 
action up to the second order in the gravitational wave tensor h^, the evolution equation 
for tensor fluctuations is obtained in Fourier space as [17] 



(A*T+ U 2 -±-h/ = 0, (3) 




where the subscript ' denotes the derivative with respect to r], the amplitude n s {rj) is 
defined by = z s h s , and the variable z s is defined by 





^.k^aWi/l-A'^, (4) 

where subscript s stands for a circularly polarized state, s =R, L. We define the right- 
handed or left-handed circular polarized state by its helicity: 

Jiyfax) = — ^ Jdk £ 4(k)/^(k)e* k - x , (6) 



=R,L 



^c^a d ef d — ke^ b , (7) 
*k c e a cd e\ d = -ke L ab , (8) 

where e f y is the polarization tensor for right-handed or left-handed circular polarization 
state. 

In equation ([3]), the important point is that the term z s "/z s depends not only 
on time, but also on the polarization mode. This readily implies that asymmetry of 
the amplitude in left- and right-handed modes may be produced, leading to a circularly 
polarized GWB. Apart from the helicity-dependent nature, the evolution equation is 
a standard form of the harmonic oscillator and one may address the quantum-mechanical 
generation of the GWB as in the case of simple inflation models. However, there exists 
a subtle issue on the break-down of linear theory arising from the singularity of the 
effective potential. Although the analysis under tractable conditions shows that the 
produced polarization-degree of the primordial fluctuations will be small [TTJ , the result 
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might not be appropriate for the practical cases. The quantitative prediction of the 
primordial spectrum may be a serious problem in the predictability of the inflation 
model. We do not discuss in details the primordial spectrum of circular polarized GWBs, 
but rather, we focus on the detectability of primordial circularly polarized GWBs. 



3. CMB power spectra from circular polarization of the GWB 

In this section, we discuss CMB anisotropies originated from a circular-polarized GWB. 
While we particularly consider the circularly polarized GWB, the linearly polarized 
GWB is shown to have no effect on the CMB power spectra due to the symmetry 
associated with spin nature of linearly polarized gravitational waves. The details are 



discussed in Appendix B 



Similar to scalar- type fluctuations, tensor-type fluctuations (i.e., GWB) cause a 
small perturbation in the photon path, producing CMB temperature and polarization 
anisotropies [28j EHl EQ]. For a temperature fluctuation map T(h) as a function of sky 
position n, let us expand it in the spherical harmonics, Yi m (n). We denote the expansion 
coefficients by aj m . Furthermore, polarization maps for the Stokes parameters Q(h) 
and U(n), which characterize the linear polarization state of the CMB, are obtained 
and are expanded by the spin- weighted harmonics Y^{n). The coefficients of these 
polarization anisotropies are decomposed into an electric part, af m , and a magnetic 
part, af m [28]. Apart from a tiny non-Gaussianity, these coefficients af m (X=T,E,B) 
are statistically characterized by Gaussian statistics with zero mean. In the case of the 
two-point statistics of CMB temperature and polarization anisotropies are completely 
specified by the six (TT, EE, BB, TE, TB ,EB) power spectra defined as the rotationally- 
invariant quantities: 



C: 



i / ( n X* n X' , „X X' 

XX' _ 1 V~V \ a lm u lm ' u lm a £m 



, \ 2 / ' [J) 



2i + 

in terms of which, 

(4l' fl L) = Cf X Su'Smm' , (10) 

where X and X' stand for T,E and B. 

Usually, the tensor perturbation produces both EE- and BB-mode polarization 
power spectra, but cross power spectra of TB- and EB-modes should vanish because 
of the parity conservation of the perturbations. However, a circularly-polarized GWB 
manifestly violates the parity symmetry, leading to non-zero values of the TB- and EB- 
mode power spectra. To be more precise, we write down the relation between the CMB 
anisotropy power spectra and the primordial power spectra of the GWB |28j : 

Cf X ' W = (4vr) 2 / edk[P iL (k) + P^(k)}AU^e(k) ; (11) 
for XX'=TT, EE, BB and TE, and 

Cj Y ' (t) = (4vr) 2 J k 2 dk[P tL (k) - P m (k)]A t Ye (k)A t Yie (k) ■ (12) 
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for YY'=TB and EB. Here, subscript w indicates the contribution from tensor mode 
and A^(/c) is photon's transfer function for X (see Appendix A). The quantities P ts (k) 
(s = L, R) are the primordial power spectra of GWB in terms of the circular polarization 
basis. The circularly polarized GWB implies P th (k) ^ P tR (k), which clearly yields non- 
zero TB- and EB-mode power spectra. 

Here, to characterize the polarization degree of GWB, we introduce the new variable 
e defined by 

P tK {k)= l -{l + e)P\k), (13) 

P th (k) = \ (1 - e) P\k) , (14) 

P\k) = P tL (k) + P tR (k) . (15) 

The variable (e + l)/2 is the fractional power of right-handed GWB with respect to that 
of total GWB. Therefore e characterizes the degree of parity violation. For instance, 
e — — 1, and 1 respectively indicate perfectly left-handed polarized, un-polarized, and 
perfectly right-handed polarized GWB. Hereafter, we simply assume that e is scale- 
independent, which might be a good approximation in the slow-roll regime |17j . 

Notice that the TT-, EE-, BB- and TE-mode power spectra remain unchanged 
irrespective of the parity violation. Thus, for CMB experiments, a measurement of the 
TB- and EB-mode power spectra is a unique probe to search for the parity violation 
in the early universe. Observationally, TB and EB-mode spectra are often used for 
a consistency null check to determine whether or not the foreground contamination is 
removed [26]. However, in our case, the non- vanishing values of the TB and EB-modes 
are essential. In this sense, detection of a circularly polarized GWB should be carefully 
investigated in practice, since the incomplete foreground subtraction may lead to a 
false detection. Nevertheless, in the next section, we will show that TB- and EB-mode 
power spectra originating from the circularly polarized GWB have some characteristic 
features, especially on large-angular scales, which might be helpful to discriminate the 
primordial origin from foreground contamination. Moreover, note that the signals of 
TB- and EB-modes power spectra originating circularly polarized GWB do not depend 
on the wavelength of CMB photon in contrast to some foreground sources. 

Finally, we comment on the TB- and EB-mode power spectra generated through 
the electromagnetic Chern-Simons term, g(x)F AF '. When the scalar field x is identified 
with the ghost or the quintessence field, this term affects the CMB polarizations after 
photon decoupling, through the rotation of the photon's polarization axis. As a result, 
we obtain non- vanishing TB- and EB-modes like Cj B = Cj E sin 2a, where a is the 
rotation angle of the polarization axis [TBI I3"T| 1321 [33] . This is even true in the absence 
of tensor fluctuations, since a non-vanishing contribution is still obtained from scalar 
type fluctuations. Thus, for a small tensor-to-scalar ratio, the shape of TB-mode power 
spectrum is essentially the same as that of the scalar-type TE-mode spectrum. In 
this respect, a non- vanishing TB-mode spectrum by the electromagnetic Chern-Simons 
term may be clearly distinguished from that produced from circularly polarized GWBs. 
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Note that the non-vanishing TB-mode is also obtained by the Faraday rotation through 
intervening magnetic fields [31]. The Faraday rotation depends on the CMB photon 
frequency [35] and it also alters the angular dependence of the TB-mode power spectrum. 
In this paper, we do not consider these two effects and just focus on the CMB power 
spectra from the circular polarized GWB. 

4. TB- and EB-mode power spectra from a circularly polarized GWB 

We now consider the amplitude and the shape of the TB- and EB-mode power spectra 
discussed in §3, taking into account the secondary anisotropies. We will show that the 
effect of reionization greatly enhances the amplitude of the TB-mode at large angular 
scales. On the other hand, the effect of weak lensing is shown to be negligibly small. In 
this and following sections, we adopt the following cosmological parameters as fiducial 
model parameters, which are taken from the best-fit values of the three year WMAP 
data (ACDM+tensor), except for the tensor-to-scalar ratio r = 0.1: 

tt h h 2 = 0.0233, tt CDM h 2 = 0.0962, fi K = 0, h = 0.787, 

r ri = 0.09, A^(0.002/Mpc) = 2.1 x 10~ 9 ,n s = 0.984, r = 0.1. (16) 

For simplicity, we assume the slow-roll consistency relation, ut = —r/8, and the 
vanishing running spectral index. The power spectra of CMB anisotropies presented 
here are all calculated based on the CAMB code [3H], with suitable modification to 
compute TB- and EB-mode spectra. 

4-1. Primary anisotropies 

Let us first show the primary anisotropies of the TB- and EB-mode power spectra 
originating from the circularly polarized GWB of a primordial origin. 

In Figure [U specifically setting the parameter e = 1 corresponding to the right- 
handed polarized GWB, the TB- and EB-mode power spectra are plotted under the 
fiducial cosmological model except for the re- ionization parameter, r r ; = 0. The results 
are then compared with the TT- and BB-mode spectra for the tensor fluctuationgEl 

Similarly to the BB-mode power spectrum, the TB- and EB-mode spectra have 
a peak at £ ~ £r, corresponding to the horizon scale at recombination (for details of 
the location of the BB-mode peak, see [36])- Also, at higher multipoles with I > 200, 
oscillatory behavior appears, which simply reflects the oscillations of the gravitational 
waves after the horizon re-entry time during the recombination epoch. A closer look 
at cross spectra reveals that while the EB-mode spectrum has many crossing points at 
higher multipoles, the TB-mode spectrum has one crossing point and the sign of the 
spectra is only changed around £ ~ 70. Further, the amplitude of the EB-mode spectrum 
is extremely small compared to the one naively expected from the BB-(EE-)mode tensor 

+ Note that the sign of the TB-mode power spectrum plotted here is opposite to the one in Ref.|13j. 
Perhaps, this differences come from the definition of polarization bases, e a ^ . Our definition follows 
that of Ref . [IT] . i.e., equations j7j and (JSj) . 
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Figure 1. The temperature and polarization cross spectra, C™^ and Cf , from 
circularly polarized gravitational waves. Here, setting the reionization optical depth, 
T r i = 0, the absolute values of the cross power spectra are plotted for fiducial model 
with e = +1. In these plots, the negative correlation is indicated by the short-dashed 
lines. As a reference, TT (dotted) and BB-modc (long-dashed) power spectra are also 
plotted. The vertical line labeled by £r(^ 100), roughly corresponds to the angular 
size of the horizon radius at recombination epoch. 



spectrum. These characteristic behaviors basically come from the projection factors in 
the photon's transfer function A^ e (k) (X=T, E, B). In Appendix C , the reasons for 
these properties are discussed in some detail. 

Apart from a tiny dependence on the density parameters such as £\ an d (e.g., 
see Ref.[37j in the case of the BB-mode spectrum), the amplitude of primary TB- and 
EB-mode spectra are mainly determined by the tensor-to-scalar ratio r and the fractional 
power of circular polarization e. In Figures [2] and El the dependence of the TB- and BE- 
mode spectra on e and r (or n^) are shown respectively. Both of the parameters e and r 
linearly alter the amplitude of spectra, but the degree of circular polarization, e, allows 
us to change the over-all sign. This is the key to discriminate the polarization states 
of the GWB. Note that in plotting Figure [31 we strictly keep the slow-roll consistency 
relation, ut = —r/8. However, the change of the spectral shape is very small and it 
would be difficult to observe it. 



4-2. Effects of secondary anisotropics 

Let us move to the discussion on the effects of secondary anisotropies generated after 
the recombination epoch. 

There are two possible sources to produce a large-angular scale anisotropy: 
reionization and the weak lensing. Among these, the weak lensing effect represents the 
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Figure 2. Dependence of circular polarization degree, e, on the TB- (left) and EB- 
modc (right) power spectra for the fiducial model except for the rcionization optical 
depth, r r i = 0. 




Figure 3. Dependence of the tensor-to-scalar ratio, r, on the TB- (left) and EB- 
mode (right) power spectra for the fiducial model with e = —1.0. In these plots, the 
reionization optical depth is set to r r i = 0, keeping the slow-roll consistency relation 
m = —r/8. 



gravitational deflection of a photon's propagation direction by the large-scale structure 
and it distorts the temperature and polarization maps of the CMB (see jl2] for a review). 
In particular, the effects from weak lensing are known as the big obstacle to detect the 
gravitational waves from the BB-mode power spectrum, since weak lensing newly creates 
the B-mode polarization anisotropy from the scalar-type perturbations, which would 
dominate over the tensor fluctuation at i ^ 100. In the case of temperature-polarization 
cross spectra, however, transformation properties of E- and B-modes do not allow the 
production of a new TB-mode correlation from the scalar-type perturbations. This is 
also true for the EB-mode spectrum. As a result, the lensing effects on the TB- and 
EB-mode spectra are negligibly small. Detailed discussion on the effects of weak lensing 



are presented in|Appendix D 
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Figure 4. Dependence of the reionization optical depth, r r ; on the TB- (left) and EB- 
modc (right) power spectra for the fiducial model with e = — 1.0. A large enhancement 
of the amplitude arises due to the re-scattering of the CMB photons during reionization. 

On the other hand, the reionization of the universe drastically changes the large- 
scale anisotropies. Although details of the reionization history are mo del- dependent and 
are currently uncertain, its effect on the CMB power spectra is mainly characterized by 
the optical depth to the beginning of reionization, r r i [401 E] • m Figure HJ we show the 
TB- and EB-mode power spectra for various values of the reionization optical depth. 
Similar to the polarization spectra of scalar-type perturbations, the resultant power 
spectra are greatly amplified and a larger value of t t \ leads to a large amplitude of TB- 
and BE-mode at lower multipoles. This is essentially the same reason as in the scalar 
TE- and EE-mode spectra that the polarization anisotropies of the CMB photon are 
newly created from a primary anisotropy by Thomson scattering at reionization. The 
power spectra are sharply peaked at large-angular scales and the peak position £ ri is 
roughly estimated as l x \ ~ ^/z^ [H]. One interesting observation is that a new zero- 
crossing point appears around I ~ 10 — 20 in the TB-mode spectrum and the amplitude 
of lower multipoles i < 6 eventually changes its sign. 

Although the precise form of the spectrum is difficult to predict analytically, the 
peak height of the spectrum caused by the reionization is roughly estimated as follows. 
First of all, the reionization reduces the fraction of photons reaching us from the 
recombination epoch. This is basically proportional to exp(— r ri ). Further, in the simple 
approximation with instantaneous reionization, the visibility function g(t]) = r' e~ T in 
equation (1A. 5[) has a sharp peak around the reionization redshift z T {, in addition to the 
primary peak around the recombination epoch. These effects explicitly appear in the 
photon transfer function A l xl or A^. Keeping this point in mind, from equation flA.lOj) . 
the transfer function for temperature fluctuation becomes 

AUk, At) ^ - / d7]e illk{v - Vo) h'e~ T « e - Tri A' NR , (17) 
Jo 

where we have only considered the dominant term. Here, A iNR represents the transfer 
function in the absence of reionization. In a similar manner, the transfer function for 



Probing polarization states of primordial gravitational waves 



11 



TB(t) 



polarization fluctuations A P , given by (IA.130 . is approximately described by 

AUvo,k,fi)= / d V e i » k ^°\- g y) * -[1 - e-^]A™ R . (18) 
Jo ID 

Here, the source function \I/ has been roughly evaluated from the monopole component 

of the temperature fluctuation as $ ~ (1/10) A^ . The prefactor, [1 — e~ Tri ], indicates 

the fractional probability of photons scattered after the reionization before reaching the 

observer, leading to a new polarization anisotropy. 

Based on the above approximations, the peak height of the power spectra is roughly 

estimated around £ ~ 2. From equations ( 1A. 1 j) ( IA79T) . we obtain 

cZ {t) - e- 2 - Q T J 2 WNR , (19) 

Cf^ -^[l-e-f^ NR , (20) 

CfT (2D 

^[l-e^}Cj^\ (22) 

where (jJ T ^ NR stands for the temperature power spectrum for tensor mode without 
reionization. For fiducial cosmological parameters, the amplitude of the TB-mode at 
£ ~ 2 is evaluated as 

cZ {t) \*4 xKT 1 | £ |(^-) [fjK 2 ] . (23) 

With an appropriate range of the reionization optical depth 0.05 < r ri ^ 0.15, the 
approximations ffT9l) - (l22l agree reasonably with numerical results of the power spectra. 

As a summary of this section, we present the full CMB power spectra, i.e., 
the combined results of the contributions from both the scalar- and tensor-type 
perturbations. Figure [5] shows the results, specifically choosing the degree of polarization 
as e = 0.1. With a slightly large value of the tensor-to-scalar ratio r = 0.1, the amplitude 
of the TB-mode spectrum becomes comparable to the EE-mode spectrum. Interestingly, 
the amplitude of the TB-mode spectrum also exceeds the BB-mode amplitude at large 
angular scales. This is even true for small degree of polarization, e ^ 0.01. 



(1 



5. Observational constraints on the circular polarization of the GWB 

Having understood the basic properties of the TB- and EB-mode power spectra, we now 
proceed to discuss the observational aspects for detecting a circularly polarized GWB. 
In section 15.11 for illustrative purposes, constraint on the degree of polarization of the 
GWB is derived based on the three year WMAP data. In practice, we must wait for 
future polarization experiments in order to get a meaningful constraint. In section I5T2"| 
future prospects for measuring the degree of polarization of the GWB are addressed 
based on a Monte Carlo analysis of parameter estimation. 
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Figure 5. CMB power spectra for the fiducial cosmology with e = +0.1, including 
both scalar and tensor mode. For comparison, three year WMAP data of TT, TE and 
EE-mode power spectra are plotted. At large-angular scale, the amplitude of the TB- 
modc (magenta, solid) exceeds the BB-modc power spectrum (green, dot-short dashed) 
and becomes comparable to the EE- mode power spectrum(6/ue ; long-dashed). 

5.1. Constraints from three year WMAP data 

Previous sections reveal that the effect of reionization largely amplifies the lower- 
multipole anisotropies and the amplitude of the TB-mode spectrum at multipoles 
£ ~ 2 — 10 would be a clear indicator to measure the polarization states of the GWB. 
While currently no definite detection of the tensor-type fluctuations has been reported, 
it is a good exercise to consider how one can constrain the circularly polarized GWB 
from lower-multipole data. For this purpose, we use the TB- and EB-mode power 
spectra taken from the three year WMAP data, currently the highest precision data at 
large-angular scales [26]. Here, particularly using the lower-multipole data of I < 16, 
we perform a global parameter estimation together with the TT-, EE- and TE-mode 
data. We use the publicly available Markov-Chain Monte Carlo code, COSMOMC [25] . 
which we modified to compute the TB- and EB-mode power spectra originating from a 
circularly polarized GWB. 

In the present analysis, we use the likelihood function for TT-, EE- and TE-mode 
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spectra available on the LAMBDA website Q. As for the TB- and EB-mode data, we 
simply assume the Gaussian likelihood function: 

£tb/eb = exp (^_^k^) , (24 ) 

with 

/£»TB/EB _ £TB/EB^ 2 

Xtb/eb = E — T 1 , (25) 

/EB 

where the quantities C and C TB / EB respectively denote the theoretical value and 
the experimental data of the TB- or EB-mode power spectra. The quantity erf denotes 
the variance of estimated power spectra at each multipole, corresponding to the diagonal 
component of the covariance matrix. Strictly speaking, the assumption (|24p is not valid 
for the three year WMAP data. For full-sky coverage, the exact likelihood function 
significantly deviates from the Gaussian likelihood function at lower multipoles [46J. 
Nevertheless, just for illustrative purpose, we adopt the Gaussian form ( J24l) . since we 
do not know the precise form of likelihood function suitable for WMAP experiment 
including TB- and EB-mode power spectra. A more rigorous treatment including the 
non-Gaussianity in the likelihood function will be discussed in the next subsection. 

To derive the constraint, we consider a spatially flat cosmology and treat the 
following eight parameters as free parameters: 

(fl h h 2 , QcDMh 2 , 9, T ri , n s , A s , r, e) (26) 

where 9 is the ratio of the sound-horizon scale to the angular diameter distance. The 
parameters ns and As are the scalar spectral index and the amplitude of the curvature 
perturbation, respectively. Then, keeping the slow-roll consistency relation, n T = — r/8, 
we perform a global estimation of the cosmological parameters. 

Figure [6] shows the constraints on the tensor-to-scalar ratio r and the circular 
polarization degree e by marginalizing over the other cosmological parameters. Top 
panel plots the projected two-dimensional contours of 68% (blue) and 95% (light-blue) 
confidence regions, while bottom panels give the marginalized one-dimensional posterior 
distribution for the parameters e (left) and r (right). Note that the constraints on the 
other cosmological parameters are also derived and our constraints reasonably match 
those obtained by the WMAP team. 

From Figure El no definite constraint on the degree of circular polarization was 
obtained. This is simply because the uncertainty in the tensor-to-scalar ratio r is still 
large. Although the 95% confidence limit of r is slightly reduced to r < 0.59 compared 
to the WMAP result with r < 0.65_$|, this is still consistent with the vanishing tensor- 
to-scalar ratio r = 0. A closer look at the posterior distribution reveals that there is a 
local maximum around e ~ — 1. However, observational errors of the TB- and EB-mode 

* http : / /lambda. gsfc .nasa.gov/product/map/dr2/likelihood_get . cf m 

jj This result is obtained using the three year WMAP data with TT-, TE- and EE-mode. Note that 
the tightest constraint is r < 0.30 with WMAP3+SDSS. 
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Figure 6. Constraints on the circularly polarized GWB from the three year WMAP 
data. Top panel shows the 68% (blue) and the 95% (light-blue) confidence regions of 
the parameters r and e. The results are obtained by marginalizing over the other 
cosmological parameters. Bottom panel shows the posterior distribution for the degree 
of polarization e (left) and the tensor-to-scalar ratio r (right). 

spectra are very large and the agreement between theory and observation is not visually 
clear. Therefore, the significance of non-vanishing e is very low. We conclude that no 
meaningful constraint on the degree of circular polarization is obtained. 

5.2. Future prospects 

Focusing on the prospects for measuring the circular polarization degree, we estimate the 
expected constraints derived from the future experiments. In what follows, assuming 
the complete subtraction of the foreground sources, we address principal aspects for 
detecting a circularly polarized GWB. We examine the two specific cases: forthcoming 
experiment by PLANCK satellite and a cosmic- variance limited experiment idealistically 
corresponding to the next-generation CMB measurement. 

As mentioned in the previous subsection, the Gaussian likelihood function for the 
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C/s is an inadequate assumption at lower multipoles [IB] and the non-Gaussianity 
arising from the cosmic variance should be properly taken into account. Further, notice 
the large cosmic variance for the TB-mode power spectrum. This is deduced from the 
diagonal component of the covariance matrix, CovJ B : 



C ° V£ - (2i+l)f sky ' (27) 



which roughly corresponds to the estimation error of the power spectrum. Here, Nj T 
and Nf B denote the experimental noises for temperature and polarization maps, and 
/sky is the fractional sky coverage. In the above expression, theoretical power spectrum 

— TT 

C e includes the contribution from both the scalar- and tensor-type perturbations. 
Thus, for a small tensor-to-scalar ratio, the dominant contribution to CovJ B always 
comes from the first term (C e + N' l F T )(C e + Nf B ) — C e C e , leading to a large 
uncertainty in the power spectrum estimation. This is true even in the absence of the 
primary TB-mode anisotropy. In this respect, definite detection of the degree of circular 
polarization requires a larger value of e and a proper treatment of the cosmic-variance 
is crucial to get the correct constraints. 

Keeping the above remarks in mind, we adopt the non-Gaussian likelihood function 



derived in Appendix E 



In " " [ JJ 

I (^TT^BB _ (Cj J' 



cTcT - (ciy 

Again, the quantities C e and C^ Y respectively denote the theoretical and the estimated 
values of the power spectra. Note that the likelihood function (128|) becomes maximum 
when C e = Cf . The above expression is the exact result for the experimental data 
with full-sky coverage / s k y = 1, but it still provides a good description for an experiment 
with almost full-sky coverage, like PLANCK. 

Based on the likelihood function ( 128|) . we perform a likelihood analysis to estimate 
the sensitivity of future experiments for constraining the parameters r and e. To do 
this, we use the TT-, BB- and TB-mode power spectra with multipoles I < 100. The 
data points for each power spectrum Cf Y are exactly set to the fiducial theoretical 
values. This is equivalent to the averaged data set over the infinite number of mock 
samples [17]. For the cosmic- variance limited experiment, we just use the form ( |28l ) 
and simply set / s k y to unity. On the other hand, for the PLANCK setup, both the 
theoretical and experimental power spectra in the likelihood function ( 1281) are replaced 
with those including noise bias contributions as cJ T//BB — > (C*J T ^ BB + A^ TT ^ BB ) and 
^yrr/BB — > (Cj T ^ BB + Nj T ^ BB ). The noise power spectra for the PLANCK experiment 
are given by 



N* x = oo^Wf 2 = (a PiX ^FWHM) 2 exp 



'beam 



(29) 
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with subscript XX being XX = TT or BB. The quantity oox is the weight factor 
per solid angle, Wg is the beam window function, and the beam size, ^beam, is given 
by 4>eam = V8 In 2/ (^fwhm) for the Gaussian beam. For the average sensitivity 
per pixel, o~p t x, and angular resolution, 9fwhm, we adopt the values for the high 
frequency instruments of 100, 143 and 217GHz channels (see Table 1.1 of (HJ for 
instrumental performance). The sky coverage of PLANCK is assumed to be / s k y = 0.65, 
corresponding to a ±20 degrees Galactic cut. 

Figure [7] displays the results for the expected sensitivity of future experiments to the 
constraint on the degree of circular polarization e for the specific tensor-to-scalar tensor 
ratio: r = 0.3 (top), 0.1 (middle) and 0.05 (bottom). In each panel, the marginalized 
68% confidence regions of the posterior distribution for e are plotted for PLANCK (red) 
and cosmic- variance limited (yellow) experiments, as a function of the true input value, 
£truc 0- At first sight, a definite detection of the degree of circular polarization seems 
difficult for small tensor-to-scalar ratios. This is simply due to the large cosmic variance 

r p r p BB _____ 

coming from the contribution C e C e (see Eqs. (l271) and (|28l)). in which the TT-mode 

— TT 

spectrum C t always gives a large value and is not much affected by the tensor-to-scalar 
ratio. From Figure d the PLANCK experiment hardly constrains the degree of circular 
polarization at r % 0.1, below which the 68% confidence level extends over the region 
fobs < and one cannot clearly discriminate between polarized and un-polarized GWBs. 

On the other hand, for the idealistic situation of cosmic- variance limited experiment, 
there still exists a window to distinguish a signature of circularly polarized GWB 
from the cosmic-variance dominated data. From Figure [7J the detectable level of the 
polarization degree can be read off: 

(», \ —0.61 
05) • (30) 

Note that this estimate is roughly consistent with the one obtained by Ref.[13], in 
which the authors reported that post-PLANCK experiment might conceivably be able 
to discriminate a value as small as e ~ 0.08 for the tensor-to-scalar ratio r = 0.7, 
comparable to our estimate of the detectable level, 0.07. However, they did not properly 
take into account the effects of reionization. Further, they only used the TB-mode 
spectrum to derive a minimum detectable e. In this respect, close agreement between 
ours and Ref.[T3] might be regarded as an accidental one. 

Anyway, a realistic value of the tensor-to-scalar ratio is expected to be much smaller 
than unity. Our results imply that a large value of e is generally required in order to 
falsify the possibility of an un-polarized GWB. However, we do not theoretically exclude 
the possibility of a perfectly polarized GWB. Though difficult, it is still worthwhile to 
explore a signature of parity violation in the universe with future CMB experiments. 

ff The upper and lower values of the 68% confidence region, [ei, £2], around the best-fit value are 
estimated from the marginalized posterior distribution P(e bs) &s J e ^ de \, s P(£obs) = 0.68 under equi- 
probability, P(£i) = P(£z)- In cases with £2 (£1) reaching 1 (—1), we simply set it to 1 (—1). 
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Figure 7. Expected sensitivity of future experiment to the constraint on the 
circular polarization degree e for the specific tensor-to-scalar ratio: r — 0.3(top), 
r = O.l(middle) and r = 0.05(bottom). In each panel, the one-dimensional marginalized 
68% confidence region of the estimated value of circular polarization degree, s bs, is 
plotted as a function of the true input value, etrue- The yellow and red shaded region 
respectively represents the expected sensitivity of the PLANCK and ideal cosmic- 
variance limited experiments. The dashed line indicates etruo = £obs- 
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6. Discussion and conclusions 

We have extensively discussed the detectability of the polarized states of primordial 
gravitational waves imprinted in the CMB anisotropies. In the early universe, the 
parity violation term originating from superstring theory or M-theory generically arises, 
which may produce a circularly polarized GWB during inflation. Such asymmetrically 
polarized gravitational waves induce a non-trivial correlation of CMB anisotropies 
between temperature and polarization modes. We have calculated the power spectra of 
CMB anisotropies generated from a circularly polarized GWB (i.e., TB- and EB-mode 
spectra). Taking into account the secondary anisotropies, we found that the effect of 
reionization creates a large amplitude of the lower multipoles of TB- and EB-mode 
spectra, which may be helpful to constrain the tensor-to-scalar amplitude ratio, r, as 
well as the degree of circular polarization of the GWB, e. We then move to discuss 
observational aspects for detecting a circular polarized GWB. Using the three year 
WMAP data, we demonstrated how one can constrain the parameters e and r from 
TB- and EB-mode data. For future prospects, we derive an expected sensitivity of 
representative experiments, i.e., PLANCK and cosmic- variance limited experiments, to 
the degree of the circular polarization. While the PLANCK experiment seems difficult 
to answer whether the GWB is polarized or not, post PLANCK experiments dominated 
by the cosmic-variance may give a meaningful constraint on the parity violation in the 
early universe. This result is interesting in the sense that the next-generation laser 
interferometers will also be sensitive to the circular polarization mode of primordial 
gravitational waves [HI 091 [50]. Although, in practice, a large value of e is required 
to falsify the possibility of an un-polarized GWB, combined results of the two different 
measurements lead to a valuable implication of the physics beyond standard inflationary 
predictions. 
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Appendix A. CMB power spectra from tensor perturbation 

In this appendix, we summarize the explicit form of the CMB power spectra for tensor 
modes. First write down the CMB power spectra as [28] : 






(X,Y =T,E,B), 



(A.l) 
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where the photon transfer functions A^(k) are the multipole moment of the function 
A t x (k,n) (see below) and their explicit expressions are given by the integral form: 



A^(fc) 



(£ + 2)! rvo 

\ \r^y.Jo M-h'e~ T + g*)PTi(x) , (A.2) 

A^(fc) = r d77(-^)^(») , (A.3) 
Jo 

KAk)= f V ° d V (-g^)P B£ (x) , (A.4) 
Jo 

with the quantity g being the visibility function defined by 

g(v) = re~ T . (A.5) 

Here, r is the optical depth for Thomson scattering between a given conformal time 
r) and the present time tjq, the quantity h is the amplitude of the gravitational waves 
and x = k(r] — 77). The prime denotes the derivative with respect to the conformal 
time 77, and the subscript t indicates the contribution from tensor modes. In the above 
expressions, the functions ^ is the source function for radiative transfer of photon and 
Pe,b£ are the projection factors for each polarization mode of photon, given by 

* = + + ^ A<T4 ~ l! Ap0 + ? Ap2 " ^ Ap4 ' (A ' 6) 

Pn(x) = , (A.7) 



.r 2 



, 2jf(a;) 4<9 x j £ (x) 



P E£ (x) = - + diUx) + -^y- + -^-L , (A.8) 

Pbi{x) = 2d xJe (x) + , (A.9) 

x 

where jt{x) is the £-th Bessel function. The expressions (1A.1|) - (1A.9I) are basically derived 
from the Boltzmann equations for photon's radiative transfer. To derive equations, first 
note that the quantities A x (k,fi) are the solutions of the Boltzmann equation, which 

are formally written as the line-of-sight integral form: 

~ , rvo 

AUk,fj)= drie~ lXfl (-ti e~ T + gV), (A.10) 
Jo 

Al(k, y) = {-12 + x\l - d 2 x ) - %xd x }A\,{k, fi), (A.ll) 
A B (k,fi) = {8x + 2x 2 dl}A t P (k,fi) (A.12) 
with the quantity A P (k, /j) being 

A P (k,n) = r drie- ix ^{-g<l>). (A.13) 



The quantities A x (k,(j,) are related to the direct observables of the temperature and 
the polarization maps, X t (h). Writing the projected maps as X t (h) = J d 3 kA x (k,n), 
the relation between A x (k,n) and A x (k,fi) are given by 

A* T (k, h) = [(1 - /i 2 )e 2 ^ R (k) + (1 - /i 2 )e" 2 ^ L (k)] A*,(fc, /i) , (A.14) 
A«.(k, n) = [(1 - /J 2 )e 2 ^ R (k) + (1 - /i 2 )e- 2 ^ L (k)] A*.(fc, /1) , (A.15) 
A* B (k, n) = [-(1 - /x 2 )e 2 ^ R (k) + (1 - /x 2 )e" 2 ^ L (k)] A|(fc, /i) . (A.16) 
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Here, the variables, £ L,R (k), are the independent random variables characterizing the 

statistical properties of the GWB. In this paper, we assume that 

(£ L *(k)£ L (k')> = S(k - k') P tL (k) , (A.17) 
(£ R *(k)£ R (k')> = S (k - k') P m (k) , (A.18) 
<£ L *(k)£ R (k')> = o . (A.19) 



Starting from the line-of-sight integral solutions of the Boltzmann equation (1A.10I) - 
(IA.12I) and using the relations flA7T4j) - flA7L9l . one can derive the expressions for the CMB 
power spectra ( lA.ll) -( lA~2l) with help of the definition and the multipole expansion 
of the anisotropies on a projected sky: 

af m = JcM YlJn) J d 3 k Ai( Vo , k, n) . (A.20) 

For details of the derivation, the readers may refer to Refs. [281 152] . 

Appendix B. Linear polarization of the GWB and CMB power spectra 

In this paper, we have mainly focused on the detectability of circularly polarized GWB. 
Here, we briefly discuss the measurability of a linearly polarized GWB. 

Let us recall that the circularly polarized states of gravitational waves are related 
to the linearly polarized states as h L = (h + + ih x )/ \/2 and h R = (h + — ih x ) / Using 
this relationship, the power spectra of linearly polarized GWB can be rewritten with 

p t+ {k) = 

= ((£ L * + £ R *)(£ L + £ R ))/2 

= P tc (k) + {P m (k) + P tL (k)}/2 

= P tc (k) + P\k)/2 , (B.l) 

= ((£ L *-£ R *)(£ L -£ R )>/2 

= -P tc (k) + {P m (k) + P tL (k)}/2 

= -P tc (k) + P t (k)/2 , (B.2) 

where we have defined the cross power spectrum between left- and right-handed 
polarized states by P tc (k) = (£ L *£ R + £ R *£ L )/2. The above expressions readily imply 
that the linearly polarized GWB comes from the non-vanishing contribution of cross 
power spectrum P tc (k). Thus, the crucial question is whether the cross power spectrum 
P tc {k) is measurable or not. 

To clarify this, consider the TT-mode power spectra. Following the definition ([9]), 
we have 

r TT(t) _ 1 y/„T* T v 

1 ~ 2£ + 1 ^ m 

= E / dQ 'J dQ J d3k 'J d " k n m (n')Y;jh) 

x (A^k'^OAT^cbM)). (B.3) 
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x 

X 



from equation (!A.20j) . The substitution of equation (1A.14j) into the above expression 
leads to 

Cj m ~E / Off dSl J d'k'J d 3 k Y em (n')Y; m (n) 

1 - /i /2 ) {e- 2i ^e R *(k') + e 2 ^ L *(k')} A^(r/ , k') 
(1 _ ft { e ^R (k) + e - 2 ^ L (k)} A* T (77o, k)) 

^yE J dn ' J dn l d3k ' J d " k Y t™(n')Y; m (n) 

x (1 - /i' 2 )(l - tf&Zto^&rfak) 

x -^ e -2i4>' + e -^ e -W) P t(k)5(k - k')/2 

+ (e- 2 ^'e- 2 ^e ii *(k , )e i (k) + e 2 ^'e 2i ^ i *(k , )e ii (k))] . (B.4) 

In the last line of equation (IB. 41) . there appears the cross-correlation of the ensemble 
between £ L and £ R , which represents the contribution from the linearly polarized GWB. 
However, further proceeding to the integral over the azimuthal angle <j), it turns out that 
this term becomes vanishing. Because of Yi m oc e imr ^, the relevant part of the integral 



It 



over can be written as 







e^e-™* = 2n5 m±2 , (B.5) 



which thus leads to 

r2n ,\ / f2n 



(7 * dfte^'e^A (7 n d^e-^A = (2vr) 2 5 m2 5 m _ 2 = . (B.6) 

Hence, linearly polarized GWB is shown to be null effect on the TT-mode power 
spectrum. Note that similar argument does hold for the other power spectra, since 
all of the photon transfer functions are written by a linear combination of e 2 *^£ R 
and e - 2i< ^ L (see Eqs.§KjA§, (lATTol) and (|ATT6|) V Thus, the cross correlation term always 
has the factor e~ 2lt ^' e~ 21 ^ or e 2j<?i 'e 2 ^, which finally vanishes after the integration over 
the azimuthal angle. 



Appendix C. Characteristic features of TB- and EB-mode power spectra 

In this appendix, we discuss the details in characteristic features of TB- and EB-mode 
power spectra originating from circularly polarized GWB. 

§4.11 reveals that while the TB-mode power spectrum has one zero-crossing point 
around £ ~ 70, the EB-mode power spectrum has many zero-crossing points with tiny 
amplitudes. These features are mainly attributed to the projection factors in photon's 
transfer function, A^ EB£ (see fjA.7[) - flA.9jl ) . 



_ ji(x) 



./•- 



P Ee (x) = -jt(x) + dln{x) + 



./•- 
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P B e(x) = 2d x j e (x) + 



4je(x) 



x 



Let us first consider the TB-mode spectrum, in which there appears the projection 
factors, Pti{x) x Pse(x), in the kernel of the integral (1A.1I) . Figure EH] shows the function, 
Pti(x) x P B i(x), as function of x = k(r] — if) for various multipoles with £ = 10 (top), 
70 (middle) and 100 (bottom). The function starts to oscillate around x ~ £ and it 
asymptotically decays as x~ 4 . Thus, the main contribution to the integral in equation 
(1A.1I) comes from the first several peaks in the oscillations and many oscillations at 
large x are almost canceled out. Just focusing on the first part around x ~ £, we 
find that the positive part of the oscillating amplitudes has relatively larger values for 
lower-multipoles, while the amplitude at higher-multipoles has slightly large negative 
amplitudes. Eventually, the values of the positive and negative amplitudes become 
comparable at the multipole £ ~ 70. These behaviors quantitatively explain the shape 
of the TB-mode power spectrum. 

Similarly, tiny amplitude of the EB-mode spectrum is also explained by the 
projection factor, Pei(x) and Pse(x). In Figure IC2j we plot the projection factor of EB- 
mode power spectrum, Pei(x) x Pbi(x) (blue dot), together with those of the BB- and 
EE-mode spectra, P-qi(x) 2 (green, short-dashed) and Pei(x) 2 (red, solid). The oscillation 
of the projection factor P B i(x) 2 is 7r/2 out of phase with corresponding one of EE-mode 
power spectrum. Thus, the amplitude of the product, Pei(x) x Pbi(x), is degraded as a 
result of phase cancellation. Hence, the amplitude of EB-mode power spectrum becomes 
much smaller than those of the EE- and BB-mode spectra. 

Appendix D. Weak lensing effect on CMB power spectra 

Here, we derive the expression for the lensed TB-mode power spectrum following the 
discussion in Ref. |I3]. The gravitational lensing effect appears as the angular excursion 
of the photon path. Since the lensing effect is only relevant at the small angular scales 
in the CMB, we consider the small scale limit. In terms of Fourier components we have 
following the expressions for the Stokes parameters: 



T(9) 



T(9 + 56) 



I 
1 



d 2 £ 



e il ^T(i) , 



Q0) 



Q(6 + 56) 



(2nf 

d 2 i 



e^ +s ^Q(£) 



(2vr) 2 




(D.l) 



where X describes the unlensed X. 

The polarization parameter Q and U can be expressed with E and B as: 



Q(£) = E(£) cos 20 £ - B(£) sin 2<f> t , 
U(£) = E(£) sin 2<f) e + B(£) cos 2(j) e , 



(D.2) 
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Figure CI. The projection factor of TB-mode spectrum, Pje(x) x Pbi(x), as function 
of x — k(r)o — r/). The top, middle, and bottom panels represent the results with 
multipolcs £ = 10, £ = 70, and £ = 100, respectively. 
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Figure C2. The projection factors of EE-, BB- and EB-mode power spectra at 
the multipole £ = 5: P^(x) 2 (red, solid), P-q^x) 2 (green, short-dash), and pEe(%) x 
PBi(x)(blue, dot). 




where (f)£ is the azimuthal angle. The ensemble average of each Fourier components 
becomes 

(X(t)Y(i)) = (27r) 2 Cp5(i- gf) (D.3) 

with X,Y = T,E and B. Using equations (ID.ll) . (ID.2[) and (ID. 31) . cross correlation 
functions C*tq and Ctu are expressed as 

Ctu ^ = /^^ cos ^^^^ (D.4) 

where we set 9 as 9 = 9 a — 9 b- The above expressions still possess the ensemble average, 
which represents the statistical average over the photon excursions caused by the lensing 
effect. In the weak lensing limit, the angular excursion is approximately described by 
the random Gaussian distribution with small dispersion. We have (e.g., Ref. [51j): 

{e iH S e A -se B)) = exp |_| [(T 2 (e) +0 os(2^)^(e)] J 

~ 1 - j[a 2 (9) + cos(2fa)a 2 2 (9)} . (D.5) 

The functions, o"g and cr|, characterize the rms fluctuations of the photon path (see 
Ref. [43 for their explicit expressions). Substituting the relation (ID. 51) into (ID. 41 ) and 
integrating over the azimuthal angle fa, we obtain 
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Ctu{9) 



2tt 



C 



TB 



Hie) |i - + {Mm + M£9)} 

(D.6) 

The above expressions finally lead to the angular power spectrum of TB mode with a 
help of the relation: 



TB 



2tt / 9d9C TV {9)J 2 (£9) 



The resultant expression becomes 



(D.7) 



TB 



(D.8) 



9d9J 2 (i9) 



-a 2 2 (9){J (l>9) 



(D.9) 



That is, the lensed power spectrum of the TB-mode is generated if and only if the 
primary TB-mode exists. No other cross spectra can create the lensed TB-mode. This 
may be explained intuitively by a simple symmetry reason; the change of TE-mode 
into TB-mode breaks parity which we do not expect from weak lensing effect. Since 
the transformation matrix Wf' is the oscillating function whose amplitude is basically 
much less than unity [3B], the lensing effect on the TB-mode spectrum can be safely 
neglected as long as the primary TB-mode spectrum is generated from the tensor-type 
fluctuations. 



Appendix E. Exact form of likelihood function 



In this Appendix, we briefly sketch the derivation of the exact likelihood function 
used in §5.21 To do this, we first follow the simplest case of the likelihood function 
with temperature anisotropy data alone. The likelihood function for the temperature 
anisotropies observed by a perfect experiment (i.e., noiseless and full-sky observation) 
has the following form: 

fr s -if 



C(T\C, 



1 



exp 



(E.l) 



where T denotes our temperature map, S is correlation matrix given by = X^(2^ + 
1)C*£ Pi{n% ■ rij)/(47r), where the Pi are the Legendre polynomials and hi is the pixel 
position on the map, and |S| denotes determinant of correlation matrix. Expanding 
the temperature map in spherical harmonics: T{n) = J2e m a ImYtm(n) , the likelihood 
function for each aj m becomes 



c(f \cT) oc n 



c, 



TT 



exp 



\n T I 2 
— TT 

2cr 



(E.2) 
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If we assume that each multipole moment aj m just follows the Gaussian statistics with 
variance of Cj T , the above expression can be reduced to ax 2 -distribution with (2£+ 1) 
degrees of freedom: 



21n£ = J2 



-(2l-l)lnC/ 1 + (21+ 1 



lnCj T 



C 



TT 



c 



TT 



(E.3) 



where Cj T denotes the estimator defined by Cj T = J2 m \ a J m \ 2 /(^^ + !)■ 

Assuming a uniform prior distribution, the posterior distribution function is 
proportional to the likelihood function as a result of Bayes' theorem. Thus, the likelihood 
function (IE. 31) can be viewed as the posterior distribution function, as a function of the 

TT 

theoretical value, C e . Then, appropriately normalizing the posterior distribution, the 
exact expression of the likelihood function for temperature anisotropy data is obtained: 

^TT* 



-21n£(cf r )=E(^+ i : 



In 



c 



c 



TT 



c7 



(E.4) 



The above result can be extended to the likelihood functions for the general case 
with temperature and polarization anisotropies. Restricting the analysis to the case of 
the temperature and B-mode polarization data, the likelihood function becomes 

a^C- 



C 



n 

em 



1 



exp 



l d 



where the vector d and the matrix C are respectively given by 
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Then, just repeating the same procedure as presented above, we obtain the likelihood 
function: 
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